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1. Introduction 

Digital images can suffer a number of different operations that take into account e.g. acqui- 
sition, coding, compression, transmission, and many others before their final use. In such 
processes, digital images might be affected by many kinds of degradations revealed by dif- 
ferent types of visual distortions. A fist of different possible types of degradations has been 



considered in the image database TID2008 by Ponomarenko et al. [!]• A key feature about 
images is their quahty. However, image quahty has not a simple definition. In a semantic 
approximation, image quality is understood as the subjective impression of how well image 
content is rendered or reproduced [2|. This definition yields to the concept of image quality 
assessment. In practice, such assessment can be achieved by two means: one is by psychophys- 
ical subjective experiments with human observers, and the other one is by objective metrics 
applied directly to digital images [s] . Image quality assessment is a very active research area 
with many contributions, techniques and models (2]. Typically, image quality can be meas- 
ured by comparison with a reference, but unfortnately such reference is not available in many 
applications. Hence, no-reference image quality assessment methods are required and they 
have been the subject of a very active research during the last years. However, a universal 
method for estimating the overall image quality is still a challenging open issue. Customarily, 
existing no-reference methods deal with simple specific degradations such as Gaussian blur 
or additive white noise. A good survey about existing no-reference sharpness metrics may 
be found in Ferzli and Karan [s] where their metric is compared to others with good results. 
Ferzli and Karan introduced the notion of just noticeable blur (JNB). It is an edge-based 
sharpness metric based on a human visual system model. Other metrics have been developed 
including blur and noise simultaneously, such as the one by Zhu and Milanfar [6], introducing 
a new concept called true image content. Their measure is correlated with noise, sharpness, 
and intensity contrast, manifested in visually salient geometric features such as edges, show- 
ing that such measure correlates well with subjective quality evaluation for both blur and 
noise distortions. However, Zhu and Milanfar technique has been designed to compare im- 
ages within the same context, while the sharpness metric of Ferzli and Karam has been 
developed to predict the relative amount of blurriness in images regardless of their context 
(note that we consider that images resulting from distorting a given original are contextual 
each other and we use the term no-contextual when the degraded images are the result of 
distorting different originals). According to Zhu, JNB technique fails to capture the trend of 
quality change in BM3D [6] denoising experiments, since it cannot handle noise well. Later 
on, Narvekar and Karam [T] proposed an improved algorithm based on the JNB paradigm 
for a no-reference objective image sharpness metric, introducing a technique they called the 
cumulative probability of blur detection (CPBD). In such work, the sharpness metric con- 
verges to a finite number of quality classes. They used the LIVE 33 database to validate 
the performance of their metric. A training-based method determines the centroids of the 
quality classes that represent the perceived quality levels. Classification is based on assigning 
the image to one of the quality classes and then using the index of the corresponding quality 
class as the metric value for that image. They include measuring experiments for Gaussian 
blur and JPEG2000-compressed images and they show that this metric performs better than 
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other known metrics. 

In summary, different reference and no-reference classifications of metrics can be consid- 
ered. According to the type of distortion they measure, specific metrics (only one type) versus 
universal metrics (all types) could be considered. According to their contextual character, 
contextual metrics, if the measures have sense only when the images have the same context 
versus no- contextual metrics, if the measure have sense by itself, regardless of the image con- 
text could be considered. In addition to the previous schemes, supervised and unsupervised 
approaches may be considered in the image quality context. 

This paper is a step forward to the approach introduced by Gabarda and Cristobal |8 
that links image directional entropy with image quality assessment. In js], the variance of 
the directional entropy of the image was introduced as a quality assessment parameter, 
presenting the highest value when the image is free of Gaussian blur or Gaussian noise and 
decreases monotonously with blur an noise increments. Gabarda and Cristobal's method may 
be classified as a no-reference, specific, contextual and unsupervised metric. In this paper we 
extend this methodology to both contextual and no-contextual applications for image quality 
assessment by introducing the von Mises distribution of the directional image entropy as the 
cornerstone of a new image quality metric for blur and noise. 

This paper is organized as follows. Section 2 gives the mathematical background required 
to understand our methodology. Section 3 presents some applications of the VMD for image 
quality measurement, along with the performances of this method against the Mean Opinion 
Score (MOS) of TID2008 image database and Differential Mean Opinion Score (DMOS) of 
LIVE image database. Results are compared with scores from other no-reference quality 
assessment existing methods. Finally, conclusions are drawn in Section 4. 

2. Description of the method 

2. A. Directional entropy 

The information content of images can be measured in different ways. One of the most 
popular is the entropy. Generally speaking, the measurement of entropy was initially proposed 
by Shannon j9j as a measure of the information content per symbol, coming from a stochastic 
information source. It can be considered as a quantifier of the uncertainty or randomness 
of a signal or, in other words, a measure of the information content of the signal. Given a 
discrete random variable X = {xi,X2, ...,xi}, the Shannon entropy of this discrete random 
variable can be computed as 



I 




(1) 



i=l 



where P{xi) is the probability of event Xj. 
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Later, Renyi 10 extended this notion to yield the generahzed entropy, whose representa- 



tion for discrete random variables is 



j=i 



Here a is a real-positive number. It can be shown that the Renyi entropy, as a generalization 
of the Shannon entropy, reverts to the Shannon entropy when a — )■ 1. 

Measures of entropy may be global, affecting the whole image, or may be localized in each 
pixel by the use of an sliding window, centered in each pixel, that gathers the neighborhood 
of that pixel with the desired geometry. 

Our proposal is to use a linear array of pixel values, centered in the interest pixel and 
oriented in a set of desired directions, in order to have a directional measure of entropy at 
pixel level. To do so, an operational sliding window L is defined, containing the + 1 pixels 
that better fit the ideal line centered in the interest pixel and directed by angle 6. N is 
required to be an even number to preserve the symmetry of the window. 

Better than considering the gray-values of the image as input to the entropy calcula- 
tion, we prefer considering the space-frequency information of the image. Space-frequency 
representations of a signal open up the possibility of exploring the image information tak- 
ing into account the position and the local frequency content of the image. One of these 



space-frequency (SF) representations is the Wigner distribution 11 . A common analytical 



framework of SF representations was introduced by Cohen 12 . The Wigner distribution 
(WD) has been selected due to its excellent properties and because it can be regarded as a 
masterform distribution function from which the other existing representation can be derived 
as filtered versions of it |13|. The WD approximation for discrete signals is better known as 
pseudo- Wigner distribution (PWD). To calculate a local PWD we will use the following 
equation 



— -1 

' 2m ^ 



W{n,k)=2 2(n + m)z*(n-m)e"^^^H^) (3) 



N 
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This approximation of the PWD is due to Claasen and Mecklembruker jl4j and is similar to 
Brenner's expression 15 . In Eq. ([s]) the variable z{n) represents the gray value of pixel n in a 



given image z. A directional PWD may be calculated by using Eq. ^ along with a pixel-by- 
pixel sliding window Dg to pick up the values z = z{—N/2), z{0), z{N/2). The central 
position of the window matches exactly pixel n on the image. Eq. ([s]) can be interpreted as 
the discrete Fourier transform (DFT) of the product r{n,m) = z{n + m)z*{n — m). Here z* 
indicates the complex-conjugate of signal z (note that z = z* for real valued images). The 
sum is limited to a spatial interval {—N/2, N/2 — 1). In Eq. (|3|, n and k represent the space 
and frequency discrete variables respectively, and m is a shifting parameter, which is also 
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discrete. By shifting the window to all possible positions over the image, the full pixelwise 
PWD of the image is produced. 

In what follows, we will take the PWD in Eq. ([s]) as the distribution to be used as input 
to Eq. ([2]). Consequently, each spatial position n will receive a certain value of entropy. 

Diverse measures of entropy may be defined upon the kind of distribution and normaliza- 
tion used with the formulation due to Renyi. Although the Renyi measures in time-frequency 
distributions formally resemble the original entropies, they do not have the same properties, 
conclusions and results derived in classical information theory. The positivity, P{n, k) > 
will not be always preserved, along with the unity energy condition X]fc-P(^'^) ~ ^^^^ ^ 
given pixel n. In order to reduce a distribution to the unity signal energy case, some kind 
of normalization should be applied. The normalization can be performed in various ways. 



leading to a variety of possible measurement definitions 16 , 17 , 18 , 19 , 20 . 

The Renyi entropy measure applied to a discrete space-frequency distribution, as the PWD 
denoted by W{n, k), can be written as 



R(n) = -hog, ^W%n,k)j 



(4) 



Here a = 3 and we have dropped the subscript a. The summation is restricted to the 
1-D neighborhood of pixel n, as described above. Again k is the frequency discrete variable. 
In order to reduce the distribution W{n, k) to the unity signal energy case, some kind of 
normalization must be done 16 . To do so, the PWD given by Eq.([3]) is normalized by 



applying the following normalizing law, as a previous step to the Renyi entropy measure: 

~ ^ Win,k)W*in,k) ^ W^n^k) 

' E{W{n,k)W*{n,k)) ZW^{n,k) ^ ^ 

k k 

Squaring of W is justified due to the real character of the PWD, and the normalizing step 
affect exclusively to index fc, when the window operation in Eq. Q is applied to pixel n; 
hence the condition ^ W{n, k) = 1 in such position is fulfilled. 

k 

It can be shown that the Renyi entropy meet the constraint < Ra{n) < logg . Hence, 
the measure can be normalized to the interval [0, 1] by applying R{n) = R{n) /log, N in Eq. 

S- 

The 1-D window used to measure entropy by this method may be oriented in any desired 
direction 9, as indicated above, and entropy is now denoted by R{n, 0). This local pixel-wise 
directional entropy can be constrained to a scalar value for the whole image by averaging all 
the R{n, 9) values by means of 

:R(0) = (i?(n,e)) = ^5^i?(n,e) (6) 

n 
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provided that M is the number of pixels in the image. 



2.B. The von Mises distribution 

The von Mises (VM) distribution belongs to probability theory and is used to handle di- 
rectional statistics for continuous probability distribution on a circular basis. It appears, in 
many respects, analogous to the normal distributions for a scalar variable. It was proposed 



by von Mises 21 to study deviations of atomic weights from integer values. This distribution 



has also been applied to diverse applications in many fields and has become an important 



tool in the statistical theory of directional data 22 



The von Mises probability density function for the angular variable 9 is given by: 

f{e\ii.^) = -^^e'^''°<'-^^ (7) 

where is the modified Bessel function of order 0, and 9 is defined in the interval [— vr, tt). 
The parameter k is responsible for how concentrated the distribution is around the mean 
direction /x. Larger values of the concentration parameter k indicate that the distribution is 
more closely grouped around the mean direction. When k = 0, this distribution is equivalent 
to the uniform distribution. This fact happens when directions in the image are endowed 
with equal probability. The uniform circular distribution is, for example, a good model for 
Gaussian noise. 

Different authors have dedicated some research to the link between images and the von 



Mises distribution. For instance, Vo et al. 24 proposed a new statistical modeling of natural 
images in the wavelet transform domain. They claimed that the von Mises distribution fits 
accurately the behaviors of relative phases in the complex directional wavelet sub-band from 
different natural images, and introduced a new image feature based on the VM model for 
image texture retrieval applications. Palacios et al. f25| presented what they called new 
tools for color image processing, based on the circularity of the hue variable of a color 
image. They give a definition of the median and the range of angular data and apply their 



results to detect hue edges. Feng 26 presented a local image orientation estimation method 
using the image gradient, using a combination of the principal component analysis (PGA) 
and the multiscale pyramid decomposition. Grana et al. |27| described a new approach to 
texture characterization for document analysis. By considering the autocorrelation matrix, 
they described image texture through a mixture of VM distributions. 

2. C. Von Mises distribution of image entropy 

As a further contribution to the applications cited in the previous section, we propose to use 
the VM distribution for modeling the directional distribution of the image entropy, calculated 
by means of Eq. The measure defined for the image by Eq. (|6| requires the use of a 
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bimodal von Mises distribution, provided that R{0) = R{6 + 7t). This requirement is fulfilled 
when we take a = 1/2, and /i2 = /^i + 7r in the following bimodal expression of the von Mises 
distribution 



/(^l/i,'^) 
from which we arrive to [231 
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2vr/o(K) 



^KCOs(9-/i2) 



f{0\fl,K) 



cosh (fi;cos(6' — fj,)) 



(9) 



27r/o(/s:) 

Eq. (|9| can be considered as the basis for modeling directional statistics of image entropy. 
One of the features to take into account in the discrete directional model for image entropy 
starting from Eq. ^ concerns the number of axis that can be defined in the image. Two 
requirements must to be fulfilled in order to consider that the directional measures, based 
on Eq. ^ are completely comparable. First all directions must be measured by windows 
consisting in the same number of pixels, and the second one is that the span of the windows 
in all possible directions must have the same length. These two constrains are only possible 
when using four directions along the axes of a regular octagon, i.e.: vr/S, Svr/S, Svr/S, Tvr/S. 
Under this requirements, we are going to restrict our calculations upon ID windows directed 
by the axes along these angles. The arrangement of the pixels in these four operational 
windows may be visualized in Fig. [T] 

The four directionalities 9i G {vr/S, Svr/S, Svr/S, Tvr/S} that we have defined are used to 
calculate four values R{6i) of entropy by means of Eq. ^ that it can be referred as vectors 
Rj = {Ri,6i). These directional entropies can be represented also by means of Cartesian 
vectors by Rj = {xi,yi)'^ = {Ri cos 9i, Ri sin 9i)^ . 

Some definitions for a resultant vector and a mean angular direction in VM distribu- 
tion have been given by Jammalamadaka and SenGupta [28] and one analysis for modeling 



circular data may be found in Bentley 29 



In our approximation we first perform a preliminary estimation of /i by a single value 
decomposition (SVD) of the directional vectors Xj, following to Feng [263, ^^^^ deals the 
image orientation estimation problem through the local image gradient. The singular vector 
Vj of the matrix 

fRicosOi i?isin6'i\ 



X 



Rf 
Hi 

T 



R2 COS 6*2 R2 sin 6*2 



-R.s sin 6'fi 



(10) 



_R3 cos 6*3 
\R4/ \i?4Cos6'4 R^sinO^/ 

corresponding to its largest singular value, gives an estimation of the parameter fi = arg(vi). 

Also, an initial estimation of k can be achieved by the following approximation due to 
Dhillon [30]: /I = 1/2(1 - R), where R = || Ri|l/4- 



7 



The second step is a gradient descent algorithm to find a better estimation for k. The 
set of values {f{6i),Ri}, for 6i = 7r/8, Svr/S, Svr/S, Ttt/S, with i = 1,2,3,4 respectively, are 
matched by a minimum square error estimation method to find the coefficients A and B 
that fit the equation Ri = Af{6i) + B. An exact fitting is attained if A = 1 and B = 0. An 
error function e{K) = \\ {A{k,), B{k,)) — (1,0) || is defined to control the k increments, in order 
to find a minimum for this error function. The algorithm runs by means of a recursive law 
Ki+i = (1 ± C) X Kj, C = 0.01 for updating the values of k. The algorithm stops when the 
error value reaches a minimum. A Matlab implementation of this algorithm may be found 
in 
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The accuracy of the estimated distribution f{6) to a true von Mises distribution is meas- 
ured by defining a fitness parameter = e~^, whose value will be 1 for the best fitting and 
for the worst. An example of estimation of the von Mises distribution for a given image 
is shown in Fig. [2] Triangles indicate the directional entropy values {Ri, 9i) used as input to 
the algorithm, along with the estimated von Mises distribution, f{9), of the image. 

To characterize the behavior of k and the algorithm has been tested with the images 
in the image database TID2008 due to Ponomarenko et al. [l]. This database contains 25 
reference color images and 1700 distorted images (25 reference images x 17 types of distortion 
X 4 levels of distortion) in bitmap format without any compression. Images are 512 x 284 
pixels in size. The test has been performed taking the 25 reference images, excluding image 
25th, provided that this last one is artificial and has been considered inappropriate for our 
study. 

Reference images have been first degraded by adding blur in iterative manner by means of 
a rotationally symmetric Gaussian lowpass filter of size 5x5 pixels with standard deviation 
a = 1, and in a second instance by adding Gaussian noise, also iteratively, with standard 
deviation a = 0.01. Considering the 24 originals and nine iterations for each original image, 
we produced a test set of 240 images. Figures [3] and |4] show the average values of n and 
if for all these images, having the degradation level as abscissa. Measurements indicated 
that (f = 0.88 ± 0.02 as the expected value of the fitness parameter for the 24 originals 
(see abscissa=0 in Fig. |4]), indicating that the fitness parameter of the VM distribution 
converges to a common value for good quality images. In general, when more blur or Gaussian 
noise is added, fitness decreases, and the variance of increases, indicating that, for severe 
degradation, the von Mises entropy paradigm for image modeling does not hold any more. 

Figures [3] to [6] show the behavior of the parameters k and (f by averaging the values from 
the outcome of our algorithm for the 24 images in TID2008 database. A hybrid blur-noise 
situation has been experimented by adding blur and noise simultaneously to the images in 
TID2008 database (see Fig. |7| and g 

Two main consequences are derived from the previous observations. First, the fitness 
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parameter has a low value {ip « 0.88) for blurred images or very noisy images, and a 
high value {ip = 0.88) for undegraded images. The second consequence is that quality of the 
images is maximum when k reaches a maximum in the series of diversely degraded images 
from a given original. 

3. Application of VM distribution to image quality assessment 

3. A. VM distribution and image quality assessment of contextual images 

The behavior of the k parameter of the von Mises distribution suggest that it can be consid- 
ered as an image quality assessmente index, when dealing with Gaussian blur or Gausssian 
noise. In order to have a quantitative evaluation of this statement, we compare the classifi- 
cation after this parameter vs. the Mean Opinion Score (MOS) for two specified distortion 
types in database TID2008. 

The image database TID2008 includes 17 different degradation types with one to four 
strength levels. The image database includes information about the MOS for each distorted 
image. The degraded images labeled as Gaussian noise, and Gaussian blur have been scored 
by means the k value determined by our VM algorithm. Later on, Kendall, Spearman and 
Pearson correlation coefficients have been evaluated by comparing our scores with the MOS 
of the images in TID2008 image database. Results are shown in Table [Tj Figures in Table [T] 
must be understood as measures in a contextual, specialized, no-reference and no-learning 
metric. 

In order to put figures in Table [T] in line with other existing no-reference methods, we 
include Table |2} where our measures are compared to the results pubhshed by Zhu and 
Milanfar [6] using the same image database and the Spearman correlation coefficient. 

Figures in Table |2] reveal that the k parameter of VM distribution has a strong correlation 
with Zhu and Milanfar measures for Gaussian blur and Gaussian noise distortions. 

3.B. Autofocusing example 

A short depth of field is an inherent limitation in optical microscopy. Basically, the focal plane 
of the microscope needs to be critically located to observe the sharper image. Hence, auto- 
focusing techniques must be included in automatized optical microscope systems. Relevant 



examples of images from bright field microscopy may be found in Valdecasas 32 . 

The behavior of the k, parameter of the von Mises distribution suggest that it can be con- 
sidered as an image quality assessmente index, when dealing with Gaussian blur or Gauss- 
sian noise. After the above quantitative evaluation of this parameter through the TID2008 
database, we present a real life example to illustrate one of the possible applications of this 
no-reference image quality assessment method for contextual images. In this example, a se- 
quence of 100 images from an optical microscope have been processed by our algorithm to 
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determine their k values. Images are 256 x 256 pixels in size. Images 10, 30, 50 and 70 of the 
stack are shown in Fig. |9} The best image according to our algorithm is number 57 in the 



sequence, that can be seen in Fig. [TOj Fig. [TT] plots the values of k for the whole sequence of 
images. 

3. C. VM distribution and image quality assessment of no-contextual images 

One of the objectives of this paper consists in determining whether the VM distribution can 
be used to construct a no-reference quality metric in the no-contextual case. The definition of 
a metric for image quality assessment based on the VM distribution for no-contextual images 
requires some kind of normalization, to assure independence of the measure from the image 
context. To do so, we have designed two possible measures taking into account the results 
in the previous section. Following these results, we postulate that both, the concentration 
parameter, k, and the fitness parameter, (p, decrease exponentially with increasing amount of 
degradation (e.g.: Gaussian blur), that can be modelled by expressions like k = noexp{—f3D) 
or ip = ipoexp{—f3D), where D is the degradation and /3 is a shape factor exclusive for each 
image. Empirically, we have found convergence at the origin for ipQ to 0.88, although there 
is not a similar convergence for ko- The second of these expressions is simpler to calculate, 
provided that hq is an unknown parameter, while po has been experimentally determined 
as ifo = 0.88 ± 0.02 (for Gaussian blur). The exponential variation of ip can be written in a 
logarithmic way as 

log(^ = logt^o - (11) 



The derivative of Eq. (11) gives d{\og(p)/dD = — /3, and /3 can be approximated by 

P = -A{\og^)/AD (12) 

for small increments, AD. 

Now, let us suppose that an arbitrary image produces an outcome parameter yjj. By 
convolving the image with an appropriate kernel, g, (i.e.: a PSF of 5 x 5 pixels and a = 1.5) 
a controlled amount of blur may be added to the image. Then, Jj+i = li * g is a little more 
degraded image than Jj. Supposing that the fitness parameter for Jj+i is v^j+i, and making 



AD = 1 for a given measuring scale, we can say, from Eq. (12), that /3 ~ \og{ipi/(pi^i) and 



introducing this value in (11) 



A = -(log¥^.-log</.o)//3 (13) 



Eq. ( 13 ) defines in this way a new no-contextual measure for image quality assessment for 



the Gaussian blur case. The measure requires the estimation of the decay constant (3 for each 



specific image as indicated in Eq. (12). This new measure will be referred here after as von 



Mises degradation measure (VMDM) and D will be expressed in VM degradation units. 
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In order to evaluate the performance of the VMDM, we have apphed our algorithm to 
the Gaussian blur images in the LIVE image database. This group is composed by 144 
degraded images plus 29 originals. For our validation test, we have selected a group of 29 
images (20 percent of the degraded subgroup), i.e.: from image number 117 to 145. Then we 
have correlated our quality scores with the DMOS scores of this group. Considering that the 
observations are positive quantities varying over many orders of magnitude, it is plausible 
to assume that the noise will be Gaussian, and the data will be well modeled as a Gaussian 
process. In this context, it is standard practice in the statistics literature to take the log of 



the data 34 . Then, we have used log(l + D^) instead of Di for correlation. Results for the 
Pearson and Spearman coefficients are shown in Table |3} 

As expected. Table [Sjpresents good peformance results, but still far from a perfect matching 
with the DMOS scores. An objective measure must necessarily have some differences when 
compared with recordings from subjective human observers. Even more, DMOS values will 
change if the database is created using different methodologies and environmental settings; 



hence, subjective scores are not comparable 35 , while the objective test will remain the same. 



However, the scores in Table [3] can be improved by applying some weighting transformation 
to the VMDM scores by incorporating to these scores some learning process that mimic 
the specific human preferences existing in the DMOS recordings accompanying the LIVE 
database. To do so, the logarithmic VMDM measures and the transformed VMDM measures 



may be related, for example, through a function / 34 given by 



D = /(Z};^) =Z} + ^aitanh(6i(D + Q)) ai,bi>0 (14) 



i=l 



where \1/ = {a, b, c}. 

Taking 1 = 5, the values for the ai,bi,Ci coefficients have been determined by a Monte 
Carlo process that maximizes the Pearson+Spearman correlation with the DMOS scores. 
We have used images 1 to 116 (the 80 percent of the 145 blurred images) in the LIVE 
database as the learning set and we have used the former set of 29 images for crossvalidating 
the measures. Our results are shown in Table HI The table includes also the results from 
Narvekar and Karam's method for the same group of distortion. 

Comparative results indicate that this method based on the parameter of the VM 
distribution is highly matching the DMOS scores for the Gaussian blur set of images in the 
LIVE database. Our results are comparable with the results given by Narvekar and Karam's 
method which is based on a sharpness metric for Gaussian blurring when applied to no- 
contextual images. Our method has the advantage to be self-contained and its performance 
can be even improved by adding a learning step. 
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4. Conclusions 

In this paper, we have introduced a new way of determining the VM distribution of the image 
information. The possible applications of the VM distribution for image quality assessment 
have been experimentally tested. The k parameter of the VM distribution has experimen- 
tally shown that can be considered as a suitable no-reference quality indicator when dealing 
with contextual images. Significative results have been presented for the Gaussian blur and 
Gaussian noise cases. Also, the defined VM distribution fitness parameter has shown its suit- 
ability as a no-reference quality indicator in the no-contextual case for blur images. At the 
same time as the relationship of the VM distribution to image quality assessment has been 
shown in this paper, this new methodology suggest new applications of the VM distribution 
of image entropy for image processing, that will be the subject of future work. 
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Fig. 1. Different configuration of pixels to define directional windows Dg, with 
N — 9 pixels and four equally spaced orientations 7r/8, Sn/S, Stt/S, Ttt/S. Note 
that periodicity is tt, provided that 9 and 9 + 7r represent the same direction. 
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Fig. 2. Von Mises distribution of the directional entropy for image 121 con- 
tained in the database TID2008 |1|. Triangles indicate the value of directional 
entropy for the target image after four different axes, n = —0.36, n = —78.75 
deg (twin pick at 101.25 deg.). Fitness parameter : ip = 0.89. 
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Fig. 3. Averaged variation of n against increasing amount of Gaussian blur for 
reference images in TID2008 database. 
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Fig. 4. Averaged variation of if (fitness) against increasing amount of Gaussian 
blur for reference images in T1D2008 database. 
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Fig. 5. Averaged variation of n against increasing amount of Gaussian noise 
for reference images in TID2008 database. 
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Fig. 6. Averaged variation of ip (fitness) against increasing amount of Gaussian 
noise for reference images in TID2008 database. 
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Fig. 7. Averaged variation of k against increasing amount of Gaussian blur and 
a fixed amount of Gaussian noise for reference images in TID2008 database. 
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Fig. 8. Averaged variation of ip (fitness) against increasing amount of Gaussian 
blur and a fixed amount of Gaussian noise for reference images in TID2008 
database. 
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Fig. 9. From left to right, capture 10, 30, 50, and 70 from a sequence of 100 
light microscope images. 
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Fig. 11. K values for the whole sequence of microscopic images. The maximum 
of K correspond to image 57 
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Tabic 1. VMDM Correlation Coefficients 



Type of noise 


Kendall 


Pearson 


Spearman 


Gaussian noise 


0.7778 


0.8052 


0.8083 


Gaussian blur 


1.0000 


0.9600 


1.0000 
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Table 2. Spearman Correlation Coefficients 



Method 


Gaussian blur 


Gaussian noise 


Blur + noise 


Zhu and Milanfar 


1.0000 


0.9760 


0.9210 


VMDM 


1.0000 


0.8083 


0.8083 
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Table 3. Correlation VMDM-DMOS in LIVE database for Gaussian blur 



Pearson 


Spearman 


0.8546 


0.9275 
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Table 4. Comparison of correlation for two methods in LIVE database for 
Gaussian blur 



Method 


Pearson 


Spearman 


Narvekar and Karam 


0.9211 


0.9449 


transformed VMDM 


0.8926 


0.9311 
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